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I. INTRODUCTION 

During the last two decades the correlation induced metal-insulator (Mott) transition has been one of the challenging 
problems in condensed matter physics' 1 . In most cases the breaking of spatial symmetry is a prerequisite for a Mott 
I , insulator 2 . The undoped high-Tc copper-oxide materials are famous example of such Mott insulator A However, the 
H ' one-dimensional Hubbard model 4 



C/3 



o 



I at half-filling represents a case where dynamical generation of a charge gap is not connected with the breaking of a 
discrete symmetry^. In Eq. we have used standard notations: n^ a — c\ a Ci_ a , U is the Hubbard on-site repulsive 
interaction, and Nij = 1 if i and j are labels for neighboring sites and equals zero otherwise. Thus the kinetic part 
' represents hops between neighboring sites and the interaction part gives contributions only from electrons on the same 
site. At half-filling the exact solution shows a uniform ground state with exponentially suppressed density correlations^ 
and gapless SU(2) symmetric spin degrees of freedomZ. This is in agreement with the large U expansion result, that 
at U ^> \t\ the model is equivalent to the effective spin 5 = 1/2 Heisenberg antiferromagnet Hamiltonian^, 

oo ; 

'HHeis = J 2, St ' + J 2, Si • Si + 2 , (2) 

o 

m ■ 

where J = At 2 /U (l — 4t 2 /J7 2 ), J' = At 4 /U 3 in the forth-order perturbation. Several very elegant mathematical tools 
have been developed for calculation of the effective spin-chain Hamiltonians in higher ordersi2*iii£. However, in 
agreement with the exact solution^, higher order terms are irrelevant and the ground state remains featureless. 

Discussions of the Mott-Hubbard transition within the framework of Hubbard type models are generally restricted 
to lattices of equivalent sites. However, the less common case where the spatial invariance of the system is broken 
I . via the introduction of two types of atoms, say "anions" and "cations", has attracted much recent interest. In 
the simplest case of a two-site ionic generalization of the Hubbard model we obtain a Hamiltonian of the two-band 
C . Hubbard model^ 

o . 

H = t^N id cl a Cj, a + —^2(n2i, a -n 2i+ i,a) 

! + Ue ^ n 2 j,t ra 2»,J. + Up ^ n 2 »+l,tn2j+l,J, ■ (3) 

_ _ I 

In the limit A = and U e = U = U Eq. © reduces to the ordinary Hubbard model. The limit U e — U = U, 
A ^ is called the Ionic Hubbard Model (IHM) 1 ^. In this paper we study the Alternating Hubbard model (AHM)A5> 
obtained from Eq. © in the limit A = 0. 

The IHM and the AHM describe two different ways how the lattice invariance can be broken. In the IHM the 
broken translational symmetry is traced already by non-interacting electrons, via the single- electron potential energy 



difference between neighboring sites (A). In marked contrast with the IHM the lattice unit in the AHM is doubled 
dynamically, via the different two-electron on-site repulsion energy on even (U e ) and odd (U ) sites. 

Current interest in the study of the Mott transition in ID models with broken translational symmetry was triggered 
by the bosonization analysis of the Ionic Hubbard Model (IHM) by Fabrizio, Nersesyan and Gogolin (FGN)i£. At U 
= the IHM is a regular band insulator with long-range ordered charge-density- wave (CDW), while in the strong- 
coupling limit, at U S> t, A it is a Mott insulator with a charge gap and gapless spin sector. The spin sector is also 
given by the same Heisenberg chain but with slightly renormalized coupling constants 1 ^ 



At 



3 i 



1 At 2 (1 + AX 2 - A 4 ) 
1-A 2 ~ IP (1 - A 2 ) 3 



4t 4 (l + 4A 2 -A 4 ) 

U 3 (1-A 2 ) 3 ' {) 



where A = A/U. 

The bosonization analysis shows that the CDW-Mott insulator transition has a complicated two-step nature^: 
With increasing U there is first an Ising type transition from a CDW band phase into a LRO dimerized phase at U£ h . 
With further increase of the Hubbard repulsion, at U^ p > U£ h a continuous Kosterlitz-Thouless transition takes place 
from the dimerized into a MI phase. The charge gap vanishes only at U — U^ h , while the spin sector is gapless for 
U > Ug p . This phase diagram was later confirmed by various numerical and analytical studiesi^iiSii^fli 21 .. Moreover, 
these numerical studies reveal a rather complex nature of the strong-coupling, U 3> A,i, insulating phase of the 
IHM. Although the spin sector of the IHM in the strong-coupling limit is qualitatively similar to that of the ordinary 
Hubbard model, the ground state is characterized by a LRO CDW pattern and therefore shows broken translational 
symmetry 18 i 19 i 20 . 

In this paper we derive the effective spin Hamiltonian in the strong on-site repulsion limit U e , U 3> t of the AHM. 
As we show below, in marked contrast with the ordinary Hubbard and the Ionic Hubbard model, in the case of 
AHM the effective spin Hamiltonian is not translational invariant and is given by a frustrated Heisenberg chain with 
alternating next-nearest-neighbor exchange 

H eff = S * ' S *+i + E( J ' - (-^OSi • S l+2 , (5) 

t i 

where 

The obtained effective Hamiltonian is the S = I shark-tooth Hamiltonian 9 ", the limiting case of which (U e = U = U) 
is a well-known result for the Hubbard model. 



II. THE STRONG-COUPLING EXPANSION APPROACH 



In this paper we apply the method developed by MacDonald, Girvin and Yoshioka in the case of the ordinary 
Hubbard chainii to obtain the effective spin Hamiltonian for the one-dimensional Alternating Hubbard model given 
by the Hamiltonian TLahm = T + V, where 

T = *E^ c L c ^ ( g ) 

i,j,ct 

V = U e ^ n 2 i,in 2 i,i + Up ^ n2t+i,t»2t+i,j. ■ (9) 

i i 

In what follows we consider the strong coupling limit, assuming U e > U 3> \t\. Contrary to the case of the ordinary 
Hubbard model where subbands can be classified by the total number of double-occupied states (doublons) Nd, in 
the case of AHM we have to deal with a system, where each band is characterized by two different numbers: the 
number of doubly occupied sites in even and odd sublattices, denoted by Nd e and Ndo respectively. The hopping term 
mixes states from these subbands. The "unmixing" of the AHM subbands can be achieved by introducing suitable 
linear combinations of the uncorrelated basic states. The S matrix for this transformation, and the transformed 
Hamiltonian, 



Heff = e lS HAHMe 



are generated by an iterative procedure, which results in an expansion in powers of the hopping integral t divided by 
on-site energies U e and/or U Q . 

This expansion is based on a separation of the kinetic part of the Hamiltonian into three terms: T\ which increases 
the number of doubly occupied sites by one, T_i which decreases the number of doubly occupied sites by one and 
To which leaves the number of doubly occupied sites unchanged. In addition, in the case of the Alternating Hubbard 
model each of these terms splits into several different terms, depending on whether corresponding hopping process 
takes place from even to odd site or vice versa. 

In particular, we split the Tq term into three separate hopping processes: 



where 



is a "hole" hopping term, 



T = T ° + T de + T do , (10) 

Tq =t Nijhi^- a cl a Cj. a hj- a , (11) 

i,j,ot 

T £ = 1 N 2ij n 2i,-c*cli, a c j,cxnj-a, (12) 



is a "pair" hopping term, which hops a pair to even site and 

T o° = t N 2i+i,j n 2i+i,- a cl l+l a c Jta nj^ a , (13) 

is a "pair" hopping term, which hops a pair to odd site. Here hi_ a = 1 — rii_ a . 

The term which increases the number of doubly occupied sites by one is also separated into two terms T\ = Tf +T°, 
where 

Tf = t ^ N 2 i, j n2i- a cl iiCt c j<a h j - a , (14) 
is the term which increase the number of doubly occupied sites by one on the sublattice of even sites and 

T l = t N 2i+l,3 n 2i+l.- a cl l+ l.a C j-ahj.- a , (15) 

is the term which increase the number of the doubly occupied sites by one on the sublattice of odd sites. 

Similarly, the term which decreases the number of doubly occupied sites by one is also separated into two terms 
T_i =Tt l +T b _ 1 , where 

T -l = 1 Y N i,2j h h-a C i,a C 2j,<x n 2j-a (16) 
i,j,ct 

and 

T°i = t Yj -^i,2j + l^i,-acJ,a C 2j'+l,Q n 2j + l,-a, (17) 
i,j,ot 

are the terms which decrease the number of doubly occupied sites by one on the sublattice of even and odd sites, 
respectively. 

One can easily check the following commutation relations: 

[y,T«] = (m + £ m , )A g T£, (18) 

where 



U q , q = e,o 

(U, 

0, q = Q 



We must emphasize that the Eq. QlSj l is true in all cases where allowed hopping processes connect only sublattices 
with different on-site repulsion. 

Let us now start to search for such a unitary transformation S, which eliminates hops between states with different 
numbers of doubly occupied sites: 

H' = e lS He- lS = H+[iS,H] + hiS, [iS,H]} + ... . (20) 

We follow the recursive schemeii which allows to determine a transformation which has the requested property to 
any desired order in t/U where U = (U e + U )/2. To proceed further we define: 

T« [{a}, {m}} = T«\T% . . . I£* . (21) 

Using Eq. (|18(l we can write 

k 

V,TW[{a},{m}]] =J2 A a i (m i + 6 mi , )T^[{a},{m}]. (22) 

i=l 

7i'( fe ) contains terms of order t k , denoted by H'w , which couple states in different subspaces. By definition [U,7i'l fc J] ^ 
and 7i'[ fe l can be expressed in the following way 

k 

{a} {m} i=l 

If in each 'k'-th order step, we choose = 5( fc_1 ) + S^ k \ where is the solution of the following equation 

[tS [k \V} = -H' [k] (24) 

and therefore equals 

5W = -< E ^ k , x r (fc) [{°>,M], 5>«/o, (25) 
„, {m} Ei=i A ^(™* + <W) i=i 

then the transformed Hamiltonian 

= e « w We -«w (26) 

contains terms of order t k /U k ~ 1 which commute with the unperturbed Hamiltonian and mix states within each 
subspace only. 



III. THE HUBBARD OPERATORS 



To treat correlations properly, it is important to know whether at the beginning or at the end of hopping process 
a particular site is doubly occupied or not. The introduction of so-called Hubbard operators^ provides us with the 
tool necessary for such a full description of the local environment. The X° fc -operator is determined on each site of 
the lattice and describes all possible transitions between the local basis states: unoccupied | 0), single occupied with 
"up"-spin I +) and "down"-spin | — ) and double occupied | 2). The original electron creation (annihilation) operators 
can be expressed in terms of the Hubbard operators in the following way: 

4,a = Xf + aXf~ a Cl , Q = X? a + aXr- 2 • (27) 

Correspondingly, in terms of creation (annihilation) operators the Hubbard operators have the form: 



xf 




X} a = -acl_ a n i>a , 








X 20 = -act c t 


(28) 


x?° 


= (1 — ^i,T)( 1 — 


X? 2 = n lA n lA 













The Hubbard operators which contain even (odd) number of electron creation and annihilation operators are Bose- 
like (Fermi-like) operators. They obey the following on-site multiplication rules Xf q Xl s — 6 qtr Xf s and commutation 
relations: 

[Xf , Xy] ± = % ( 5 qr Xf ± S ps X^), (29) 

where the upper sign corresponds to the case when both operators are Fermi-like, otherwise the lower sign should be 
adopted. 

It is straightforward to obtain that 

T o° = ( EE^^ 0l ° Q ' ( 3 °) 

i,j a 

T o° = t E E N^Xl^xp 2 , Tt = *EE N^X'rx-* 2 (31) 

i 7 j a i,j a. 

T ? = 1 E E aN 2i+ltj xi7 + «x°<*, zy = * * N *J x ir a Xj a (^) 

Tli = ^EE "Nw^XfXg* Tt 1 = t ^2 E oiN^iXfX^ (33) 

i.j a i,j a 

One can easily find that the spin S = 1/2 operators can be rewritten in terms of the X-operators in the following 
way 

St = 4,^,1 =X+~, Sr = c+^-.t = Xr+ , S f = \{Xf + - X~ ) . (34) 

IV. EFFECTIVE HAMILTONIAN IN THE HALF-FILLED BAND CASE 

In what follows we focus on the case of the half-filled band. In this particular case the minimum of the interacting 
energy is reached in the subspace with one electron per each site. Therefore, no hops are possible without increasing 
the number of doubly occupied sites and for any state in this subspace | ^ls) 

Tt x | = T°_ x | * LS ) = To | * LS ) = . (35) 

Equation Ij35(l may be generalized to higher orders 

T k [m]\^ LS ) = 0, (36) 

if 

k 

M^[m]=J2 m i<° ( 37 ) 

i—n 

for at least one value of n. Equation (|36J) can be used to eliminate many terms from the expansion for 7i' in the 
minimum (V) subspace. Thus, in the fourth order of T, the perturbed Hamiltonian has the form: 

H ,(i) = _ 1 T o iT o _ 1 ye^e 
U o U e 

j" rpQ rpQrpQrpO ' rpe rpderpQrpO ^ rpe rpQrpde.rpO ^ rpQ rpdorpderpO 

_ _ rpe rpOrpQrpe ^ rpQ rpdorpOrpe ^ rpQ rpQrpdorpe ^ rpe rpderpdorpe 

- -4r , T° , T?T? + X;T° T°T° ,T?+( — \ — + — \ — ) T° , T°T e , Tf 
+ ^T^TZT^TZ - ^T^TUnn + + ^-)r« 1 T 1 e T° 1 T 1 ° 

^ rpO rpe rpOrpe rpe rpO rpQrpe 



u u e (u + u e ) -i-iii u*(u„ + u a y 

' -T^T^TfT? - 1 T° 1 Tl 1 TiT° (38) 



U U e (U + U e ) -i-iii U*(U + U a ) 



Using Eqs. (|30|) - (|34[1 . one can easily rewrite the products of T-terms in (|38|) via the Hubbard X operators. Using 
the spin 5=1/2 operators in Eq. (|34[1 one obtains: 

T O iT o = T e iT e = _ 2t 2 ^ ( g. . ^ _ (39) 

i 

T° x TqTqTi = T^T^T^T? = T1 x TqTqTi = Tt x T^T^ Tl = -2i 4 ^(S* • S i+ i - |) (40) 

i 

TliT de T T o = T O iT T do T e = ^ (g . . g. + i _ ^ + ^ ^(S 2l+1 ■ S 2i+3 - J) (41) 

i i 

TliT T de T o = T O iT do T T e = _ 2j A ^ ( g. . ^ _ lj + 2t ^(S W • S 2l+2 - 1) (42) 
T o_ iT o_ iT o T o = T e iT e iT e T e = g f 4 J- (Si ■ S i+1 - ^) (Sj • Sj+i - \) (43) 

= 4i 4 ^(S, • S i+1 - i)(S, ■ S j+1 - \) (44) 
= At 4 J2 (Si-Si+i-^S^S^-^) (45) 

l,z,i+l 

Therefore we find that in the 4 th order approximation the strong-coupling effective spin Hamiltonian for the Alter- 
nating Hubbard model is given by 

We// = ■ S »+l ~ 7) + J l E( S 2'+1 ' S 2'+3 - 7) + J2 ' S ^+2 ~ \), (46) 



where 
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= 2i 4 (A- 
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u -u e 


h t/f 





4*)> (47) 



(48) 
(49) 

The effective Hamiltonian thus obtained is that of a frustrated Heisenberg chain with alternating next-nearest-neighbor 
exchange — . Note that the nnn exchange is larger for two spins separated by a site with low on-site repulsion than for 
spins separated by a site with high on-site repulsion. 



V. CONCLUSION 



In this paper we have derived the effective spin Hamiltonian which describes the low-energy sector of the one- 
dimensional half-filled alternating Hubbard model in the limit of strong on-site repulsion. The effective spin model is 
given by the Hamiltonian of the Heisenberg chain with alternating next-nearest-neighbor exchange. This model was 
intensively studied in the last few yearsSi2iS5iSSiSl. Unfortunately, conflicting results have been reported in these 
studies regarding the relevance of the alternating nnn exchange. In some studies^iSSiSl it was claimed that in the 
limit of small frustration J'/ J < 0.5, the asymmetry of the nnn exchange destabilizes the isotropic Heisenberg fixed 
point and leads to a new phase with gapless excitation spectrum and vanishing spin-wave velocity. However, other 
studies^*2& claim that the alternating nnn exchange is an irrelevant perturbation. We believe that detailed numerical 
studies of the phase diagram of the Alternating Hubbard model may shed more light on this topic. The weak-coupling 
bosonization analysis of the same model is in progress and will be published separately. 



VI. ACKNOWLEDGMENTS 



It is our pleasure to thank Henrik Johannesson, Arno Kampf, Michael Sekania and Irakli Titvinidze for many 
interesting discussions. GIJ also thanks Dionys Baeriswyl for kind hospitality and many interesting discussions 
during his stay at the University of Fribourg, where part of this work has been done. We also acknowledge support 
by the SCOPES grant N 7GEP J62379. 



Appendix 



Using the technique developed above it is straightforward to derive the strong coupling effective Hamiltonian in the 
case of the Alternating Ionic Hubbard model, Eq. ©, in the limit A ^ and U e ^ U ^> A, t. 

In this case the effective Hamiltonian is also given by the Heisenberg Hamiltonian with alternating next-nearest- 
neighbor exchange, Eq. (|46[) . but with the following exchange coupling constants: 



J = 2r 



2f 



1 



U - A U e 
4 



A 

4 



1 



((7 e + A)3 ((7 D -A)3 U e (U e + A) 2 (U ~A) 2 U 
4(C/ e + U ) 1 1 



(U e + A) 2 (U a - A) 2 (U e + A) 2 U U e (U D -A) 2 
2 2 



Ji 



(A + U e )(U D - A)U Q U e (A + U e )(U - A). 
1 , 1 U e + U 



2t 



(C/ e + A)3 ((7 D -A)3 (U e + A) 2 (U - A) 2 U (U e + A){U - A) 
1 1 U e + U 2 



(U e + A) 3 (C7 -A) 3 (U e + A) 2 (C7 - A) 2 U e (U e + A)(U a - A) 



(50) 
(51) 
(52) 



Note, that in the case when the ionic term is large on the site with larger on-site repulsion, i.e. U e > U a and A > 0, 
the asymmetry in the nnn exchange in the case of AIHM 



SJ'(A ^ 0) = 2t 4 



(U e ~ U ) 



U e U (U e + A)(U a - A) 



is larger than it is in the case of AHM 



SJ'(A = 0) = 2t 



A ( Ue-Uo) 



(53) 



(54) 



for A > (A < U a ,U e ). 
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